Bound states of neutral particles in external electric fields 
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Neutral fermions of spin ~ with magnetic moment can interact with 
electromagnetic fields through nonminimal coupling. The Dirac-Pauli 
equation for such a fermion coupled to a spherically symmetric or cen- 
tral electric field can be reduced to two simultaneous ordinary differential 
equations by separation of variables in spherical coordinates. For a wide 
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values can be found analytically. The degeneracy of these energy levels 
turns out to be numerably infinite. This reveals the possibility of con- 
densing infinitely many fermions into a single energy level. For radially 
constant and radially linear electric fields, the system of ordinary differen- 



^ ■ tial equations can be completely solved, and all bound-state solutions are 

obtained in closed forms. The radially constant field supports scattering 
solutions as well. For radially linear fields, more energy levels (in addition 
to the critical one) are infinitely degenerate. The simultaneous presence 
of central magnetic and electric fields is discussed. 
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I. INTRODUCTION 



In relativistic quantum theory, a charged fermion of spin \ moving in a background electro- 
magnetic field is described by the Dirac equation with minimal coupling to the vector potential. 
In 1941, Pauli extended this equation to include an additional nonminimal coupling term which 
takes into account the interaction caused by the anomalous magnetic moment of the charged 
particle [1]. This extended equation is usually called the Dirac-Pauli equation. Many works 
have been devoted to the investigation of exact solutions of this equation in various electro- 
magnetic fields, say, a constant uniform magnetic field, an electromagnetic plane wave, and 
more complicated ones [2]. A constant central (spherically symmetric) electric field was also 
considered by some authors. In this case the Dirac-Pauli equation is separable in spherical 
coordinates, however, exact solutions of the radial equation have not been found in closed 
forms [2]. In these studies, the nonminimal coupling is conceptually taken as some correction 
to the minimal one (though the correction is considerable for protons), and the simultaneous 
presence of both couplings causes some mathematical difficulty 

In this paper we consider neutral fermions of spin | with magnetic moment. Without 
electric charges, such particles can still interact with electromagnetic fields through nonminimal 
coupling, and can be well described by the Dirac-Pauli equation. On the one hand, without 
the minimal coupling, the Dirac-Pauli equation is simpler. On the other hand, the interaction 
solely from nonminimal coupling has not attracted enough attention, especially before the 
discovery of the Aharonov-Casher (AC) effect [3-5]. Since the AC effect is a consequence of 
the nonminimal coupling and has been observed in experiment [5] , one may become interested 
in other consequences of the interaction. For instance, it may be of interest to study bound 
states of neutral fermions in external electromagnetic fields, especially when exact solutions 
are available. It appears that this problem was not considered in the literature as far as we 
know. The purpose of this paper is to deal with this problem. It is organized as follows. 

In the next section we consider the Dirac-Pauli equation of a neutral fermion of spin |, with 
mass m n and magnetic moment fi n , interacting with an external electromagnetic field through 
nonminimal coupling. For spherically symmetric or central electromagnetic fields, it can be 
shown that the total angular momentum is a constant of motion. By separation of variables 
in spherical coordinates, the stationary Dirac-Pauli equation in a central electric field, which 
involves four partial differential equations, can be reduced to a system of two coupled ordinary 
differential equations (ODE) for two radial wave functions. Given a specific electric field, one 
can in principle solve the system of ODE to obtain the radial wave functions and determine the 
energy levels for bound states. For a wide variety of electric fields, one can find bound-state 
solutions of critical energy value m n or — m n in analytic forms. It turns out that these critical 
energy levels are infinitely degenerate. This is interesting because it reveals the possibility of 
condensing infinitely many fermions, say, neutrons, into a single energy level. Electric fields 
that support a finite number of critical bound states are also discussed. In Sec. Ill we study 
a radially constant field, in this case the system of ODE can be completely solved, and we 
have scattering as well as bound-state solutions. All bound-state solutions are given in closed 
forms. Only the critical energy level has infinite degeneracy. In Sec. IV we deal with radially 
linear electric fields. The system of ODE is also completely solvable. In this case we have 
only bound-state solutions, and many of the energy levels are infinitely degenerate. In Sec. 
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V we discuss the simultaneous presence of central magnetic and electric fields. In this case 
separation of variables is still possible in spherical coordinates. But the reduced system of 
ODE involves four coupled equations for four radial wave functions, and is thus much more 
difficult to solve. Some other remarks and discussins, say, the nonrelativistic limit, are also 
included in this section. 



II. NEUTRAL FERMIONS IN CENTRAL ELECTRIC FIELDS 

We work in (3+l)-dimensional space-time and use the natural units where h — c — 1. 
Consider a neutral fermion of spin \ with mass m n and magnetic monment /i n , moving in an 
external electromagnetic field described by the field strength F^ u . The fermion is described by 
a four-component spinorial wave function \1/ obeying the Dirac-Pauli equation [2, 6] 

{i-fd^ - \n^ v F^ - m n )tf = 0, (1) 

where 7 M = (7°,7) are Dirac matrices satisfying 

{-frf} = 2gT (2) 

with = diag(l, —1, —1, —1), and 

^ v = \b^n (3) 

It can be shown that 

\o^F^ v = ia ■ E - £ • B (4) 

where E is the external electric field and B the magnetic one, a = 7°7, and S fe = ^e ktj a lj 
where e faj is totally antisymmetric in its indices and e 123 = 1. If both E and B are independent 
of the time t, one may set 

*(t,r) = e - i£ V(r), (5) 
and obtain a stationary Dirac-Pauli equation for ip: 

= S^, (6a) 

where the Halmiltonian H is given by 

H = a ■ p + ijvy • E - /i„7°S • B + 7°m n , (6b) 

where p = — iV or p k = —idk- 

Now let us consider spherically symmetric or central fields 

E = E(r)e r , B = B(r)e r , (7) 
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where r is one of the spherical coordinates (r, 9, <fi) and e r is the unit vector in the radial 
direction. As usual we define the the orbital angular momentum L = r x p. It is not difficult 
to calculate the commutator [L, H] and it turns out that [L, H] ^ even for free particles. 
For central fields, however, it can be shown that the total angular momentum J = L + S 
where S = |S is a conserved quantity, i.e., [J,H] = 0. Thus one can have a common set 
of eigenstates for (H,J 2 ,J Z ). Because S 2 = | is a constant operator, it is also a conserved 
quantity. Unfortunately, L 2 is not conserved and cannot have a common set of eigenstates 
with (H, J 2 , J z ). If B(r) = 0, we have a further conserved quantity K = 7°(S • L + 1) which 
commutes with both H and J. In this case one can have a common set of eigenstates for 
(H,J 2 ,J Z ,K,S 2 ). 

To solve the Dirac-Pauli equation, one should choose a specific representation for the 7 
matrices. Here we use the Dirac representation [6]. In this representation we have S = 
diag(cr, er) where cr are Pauli matrices, and J = diag(j, j) where j = L + |<x is a 2 x 2 matrix. 
In this section we only consider a central electric field. The simultaneous presence of a central 
magnetic field will be discussed in Sec. V. We define ip = (</?, x) T , where r denotes transpose, 
and both and x are two-component spinors. The stationary Dirac-Pauli equation (6) then 
takes the form 

<j ■ (p - i/i n Ee r )if = (£ + m Q )x, (8a) 

cr • (p + i/j, n Ee r )x = {£ - m n )ip. (8b) 

Here four partial differential equations are involved. We are going to simplify these equations 
by separation of variables in spherical coordinates. Let us define the two-component spinors 



( /Z + m + l \ 
-Y lm (0,<p) 




21 + 1 



m 



21 + 1 



Y, 



l,m+l 



1 = 0,1,2,... ;m = -(/ + !),-/,...,/; 



(9a) 



( /Z-m+1 \ 
Yi+i, m {9, (p) 



21 + 3 
ll + m + 2 



Y, 



l+l,m+l 



(OA) 



, 1 = 0, 1,2,... ;m = -(/ + !), -/,...,/. (9b) 



V V 2/ + 3 

Here Yi m (9,<j>) are spherical harmonics as defined in Ref. [7]. Both of them are common 
eigenstates of (j 2 , j z , L 2 , S 2 ) with eigenvalues 



((Z + i)(/ + §), m+i,Z(Z + l),f) 



and 



((j + !)(j + §), ro +!,(j + i)(j + 2),f), 

respectively. It can be shown that 



lm\ 



(10) 



and 



^fL(o,<P) = ifL(OA), 



(lla) 
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<r-Lfr m (e, ( f ) ) = -(i + 2)fr rn (e, ( f ) ). (lib) 

The relation 

i 

a ■ p = — i(a ■ e r )d r + -(a ■ e r )(cr • L) (12) 

is also useful in the following. With these preparations we can simplify Eq. (8) for two different 
kinds of solutions. 

The first kind of solution to Eq. (8) is ip + = (y? + , x + Y \ where 

y? + (r, 9, 0) = « + (r)/+ (9, 0), X +(r, 9, 0) = iv + (r)ff m (9, 0). (13) 

Note that -0 + is a common eigenstate of (J 2 , J Z ,K, S 2 ) with eigenvalues ((/ + \)(l + §),m + 
I + 1, |), but it is not an eigenstate of L 2 . Using the relations (10-12), it is not difficult to 
show that Eq. (8) now reduces to a system of first-order ODE for the radial wave functions 
u + (r) and v + (r): 

du + I 

— h fi n Eu + u + = -(£ + m Q )v + , (14a) 

dr r 

dv + I + 2 

— a n Ev + + v + = (£- m n )u + . (14b) 

dr r 

Because 9 and are not defined at the origin, the appropriate boundary conditions for u + and 

v + are 

|m+(0)|<oo (1 = 0), m + (0)=0 (Z^O), (15a) 

V +(Q) = V/. (15b) 

Of course they should also satisfy appropriate boundary conditions at infinity. For bound-state 
solutions to be considered below, they should fall off rapidly enough when r — > oo such that 
ip + is square integrable. For scattering problem, they should be finite at infinity. Given a 
specific form for E(r), one can solve Eq. (14) at least numerically. This is much simpler than 
dealing with Eq. (8). For £ ^ — m Q , one can express v + in terms of u + by using Eq. (14a), 
and substitute it into Eq. (14b) to obtain a second-order ODE solely for u + : 



d 2 u+ 2 du + 
dr 2 r dr 



£ 2 -< + ^ - H 2 n E 2 + 2(1 + IK- - /(/ + 1} 



dr 11 r r 
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u + = 0. (16) 



This is similar to the radial Schrodinger equation in a central potential. It can be exactly 
solved for some specific form of E(r). This will be studied in the subsequent sections. When 
Eq. (16) is solved, it is easy to obtain v + . If £ = —m n , one can directly solve Eq. (14) without 
difficulty. 

The second kind of solution to Eq. (8) is vp~ = (tp~, X~Y ■, where 

cp-(r,9, ( j ) )=u-(r)fr m (9, ( f>), X "(r, 6, 0) = iv~(r)f^ (9, 0). (17) 

Note that -0~ is also a common eigenstate of (J 2 , J z , K, S 2 ) with eigenvalues ((/ + 1)(/ + 1), m + 
|, — (I + 1), |), As before, Eq. (8) reduces to a system of first-order ODE for the radial wave 
functions u~(r) and v~(r): 

^ — h fi n Eu~ + ^—^-u~ = -(£ + m n )v ", (18a) 
dr ' r 
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— u n Ev v = (£ — m n )u . (18b) 

dr r 

This is similar to Eq. (14). If £ ^ m n , one can solve Eq. (18b) for u~, and substitute it into 
Eq. (18a) to obtain a second-order ODE solely for v ~: 



d 2 v- 
dr 2 



+ 



2dv 



r dr 



+ 



^ - ^E 2 - 2(1 + l)/i n - - 

ar r r z 



0. 



(19) 



This is similar to Eq. (16). Note that the appropriate boundary conditions for u and v at 
the origin are 

u"(0) = V/, (20a) 

|u _ (0)|<oo (1 = 0), u"(0) = (Z^O). (20b) 

Thus Eqs. (16) and (19) have the same boundary conditions at the origin. Also note that 
they interchange if E(r) — > —E(r). If £ — m n , Eq. (19) is invalid, and one can solve Eq. (18) 
directly. 

Using the completeness relation of the spherical harmonics, it can be shown that the two- 
component spinors f^(0,(j>) and f lm (0,(fi) constitute a complete set on the sphere. More 
specifically, we have 

oo I 

E E [^(^0)/i t (^0O+/^0)//™ t (^0O]=%os^-cosW0-0 / )- (21) 

1=0 m=-(/+l) 

Therefore all possible forms of solutions to Eq. (8) are included in Eqs. (13) and (17). 

In the subsequent sections we are going to solve Eqs. (16) and (19) for radially constant 
and radially linear electric fields. Before dealing with these specific cases, we would like to give 
some special bound-state solutions for more general forms of the electric field. We assume that 
E(r) behaves like r~ 1+Sl when r — > and like r~ 1+&2 when r — > oo where 8 1 and 5 2 are positive 
numbers, and is regular everywhere except possibly at r = 0. If fi n E(r) > for r > r + where 
r + is some finite radius, we have the following solution to Eq. (14) with energy level £ = m n . 



u, 



r rr 

A^r'exp - / fi n E(r')dr' 
L jo 



0, 



(22) 



where is a normalization constant. This obviously satisfies the boundary conditions (15) 
at the origin. It is a bound-state solution because it falls off rapidly enough to be square 
integrable. It is remarkable that the energy eigenvalue does not depend on the quantum 
numbers / and m (or j — I + | and rrij = m + |). Thus the degeneracy of this energy level is 
numerably infinite. This is somewhat similar to the situation of a charged particle in a magnetic 
field with infinite flux in two dimensions [8]. To our knowledge other similar situations were 
not encountered previously in the realistic three-dimensional space. This is interesting because 
it reveals the possibility of condensing infinitely many fermions, say, neutrons, into a single 
energy level. If fj, n E(r) < for r > r_ where r_ is some finite radius, then we have the 
following solution to Eq. (18) with energy level £ = —m n . 



u, 



"(r) = 0, v l (r) = A l r'exp 



H n E(r')dr' 



(23) 
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where Aj is a normalization constant. This is also a bound-state solution, and the energy 
level is infinitely degenerate. Here we have a negative-energy solution, and will have more in 
the following sections. The presence of negative-energy eigenvalues and eigenstates is a quite 
general feature of relativistic quantum mechanics. Though these solutions are nonphysical in 
the one-particle theory, it is well known that they correspond to antiparticles after second 
quantization. Therefore we will not exclude these solutions in this paper. 

If E(r) ~ k/t for large r where k is a constant, the situation is of special interest. In 
this case the nonvanishing component of the critical solutions [uf(r) for ji n K > or vf(r) for 
fi n K < 0] does not fall off exponentially at large r, but behaves like r l ~^ nK \. To be normalizable, 
one should have / < |/i n /t| — §. Therefore to have at least one critical bound state, one should 
have \fJ, n K,\ > |. When \/i n K,\ — | is a natural number, we have \/i n K,\ — § critical boumd states 
(degeneracy over m has not been taken into account). When | not a natural number, 

the number of critical bound states is [|/i n /t| — |], where the square bracket denotes the integral 
part of a number. This is similar to the situation of a charged particle in a magnetic field with 
finite flux in two dimensions [8]. 

For a specific electric field, critical bound states with £ = m n and those with £ = —m n 
do not appear simultaneously. This spectral asymmetry is also similar to that of a charged 
particle in a magnetic field in two dimensions. Therefore vacuum polarization similar to those 
in two dimensions for charged particles [9-11] or neutral ones [12] may be expected for the 
present system after second quantization. 

To conclude this section we write down the normalization condition for bound-state (or so 
called square integrable) solutions: 

dr^ ±t (r)^ ± (r) = jf [u ± (r)] V dr + J [^(r)] V dr = 1. (24) 

The normalization constants in Eqs. (22) and (23), and those in the following sections are to 
be determined by this condition. 



III. RADIALLY CONSTANT ELECTRIC FIELDS 



In this section we consider radially constant electric fields E(r) = E where E is a constant. 
As pointed out before, Eqs. (16) and (19) interchange when E(r) — > —E(r). So we need only 
consider a positive E or a negative E . For convenience we assume that /J, n E > 0. Now Eq. 
(16) takes the form 



d 2 u + 2 du^ 
+ 



dr 2 r dr 



+ 



£' 



2f2 , 2(l + iy n £ 1(1 
A^nA) H 



1) 



u 



= 0. 



(25) 



This has the same form as the radial Schrodinger equation in an attractive Coulomb field. 
The difference is that the "Coulomb field" (the next to the last term in the square bracket) 
here depends on the quantum number I. Thus the energy levels will depend on I as well as a 
principal quantum number or a radial quantum number. As Eq. (25) is familiar in quantum 
mechanics, we will give the solutions only. Remember that Eq. (25) is invalid for £ = — m n . 
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We have the bound-state energy levels 



£01 = ™> n , n r = 0, 



(26a) 



£ 



n r l± 



± 



• 2 (n r + i + iy-(i + iy 
m » + ^° ' (n r + i + iy 



n r 



1,2 



(26b) 



Here n r is a radial quantum number. When n r = we have a positive critical energy level 
given in Eq. (26a). Though it is independent of I, we keep the subscript I to make a clear 
correspondence to the corresponding wave functions below. For n r ^ we have positive- and 
negative-energy levels, indicated by the subscript ± in Eq. (26b). The corresponding radial 
wave functions are 

ul rl± (r) = A nrl Je-^Ll+\p\ (27a) 

A'n-E'o 



V tl±( r ) = A nrl±- 



for n r 7^ 0, and 



n r + I + l)(£ nrl ± + m n ) 

«J(r) = V^, 
4(r) = 



^ +1 e^ 2 e 3 i(p) 



for n r = 0, where 



2(1 + l)fi n E 
(n r + I + 1) ' 



(27b) 

(27c) 
(27d) 

(28) 



and L^ r +1 (p), etc., are Laguerre polynomials defined in Ref. [13], which are different from those 
used in Ref. [7]. Note that the superscript + indicates the first kind of solutions (13), while 
the subscript ± indicates the sign of the energy levels. It is seen from Eq. (27) that negative- 
and positive-energy solutions have the same functional form, but the coefficients are different. 
The normalization constants are 



[n r + i + iy 



2(1 + l) 3 n r ! 



[n r + 2l + 1)! 



n r l± 



£ 



n r /± 



for n r 7^ and 



A i = 



(29a) 



(29b) 



for n r = 0. The degeneracy of the energy level £ nr i+ or £ nr i- is 21 + 2. As the energy level 
£qi = m n is actually independent of /, its degeneracy is numerably infinite. Indeed, the solution 
(28) is a specific realization of the solution (22) discussed before. 

When £ = — m n , Eq. (25) is invalid. In this case one should deal with Eq. (14) directly. 
It is easy to show that this energy value corresponds to a trivial solution. Thus all first-kind 
solutions are included in Eq. (27), and the corresponding energy levels are given by Eq. (26). 
Note that all energy levels have absolute values less than \jm\ + fi^E^. When \£\ exceeds this 
value, we have scattering solutions to Eq. (25). This will not be discussed here. 

Now we turn to Eq. (19), which in the present case becomes 
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d 2 V~ 



2dV 



dr 2 r dr 



+ 



'2 2 2z?2 2(l + l)fi n E 1(1 + 1) 
£ - m n - fi n E — — 



0. 



(30) 



Since /i n E > 0, this is equivalent to the radial Schrodinger equation in a repulsive Coulomb 
field. In this case only scattering solutions are available. These scattering solutions have energy 

Eq. (30) is invalid. Then we may deal 



£ > ^Jm 2 n + ii\El or E < -\jm 2 n + ulEfi. If £ = m n 
with Eq. (18) directly. It turns out that this energy value corresponds to a trivial solution. 
We thus conclude that there is no bound state of the second kind in the present case. 

To finish this section we estimate the "Bohr radius" of the neutron in this radially constant 
field. It is roughly equal to a~^. For the critical-energy state, n r = 0, and a^ 1 = (2fi n E )~ 
In the MKS system it reads 



a. 



-1 
0/ 



he 2 



2fi n E 



(31) 



We take |£o| = 5.15 x 10 11 V/m, the electric field strength at the Bohr radius of the hydrogen 
atom. For neutrons, we have a^ 1 = 9.5 x 10~ 4 m. This is a macroscopic length scale but 
rather small. However, it might be not easy to realize a radially constant central electric field 
with the above magnitude in the laboratory. We do not know whether there exists some such 
field somewhere in the universe. 



IV. RADIALLY LINEAR ELECTRIC FIELDS 

In this section we turn to another exactly solvable field, the radially linear electric field 
E(r) = f3r where f3 is a constant. The electric charge density that produces this field is 
p c = 3(3 /4n in the Gaussian units, which is a constant. To realize the above central field, 
however, the electric charge density should become zero outside some large sphere where the 
particle under consideration cannot reach practically. Otherwise the electric field would be 
zero everywhere. In the region of interest (inside the large sphere) the field is then radially 
linear. For reasons given earlier, we need only consider one sign of (5. For convenience we 
assume that f3pL n > 0. Eq. (16) then takes the form 



d 2 u + | 2du+ 
dr 2 r dr 



£ 2 -m 2 n + (21 + 3)/^ n - (3 2 fi 2 n r 2 - 



1(1 + 1) 



u + = 0. (32) 



This is not valid for £ = —m Q . In the latter case one can solve Eq. (14) directly and obtain a 
trivial solution. Thus all nontrivial solutions of the first kind are those arise from Eq. (32). The 
equation (32) has the same form as the radial Schrodinger equation for an isotropic harmonic 
oscillator. The difference is that the "energy" here depends on the quantum number /. Thus 
the dependence of the energy levels on the quantum numbers will be different from that of the 
isotropic harmonic oscillator. Since Eq. (32) is also familiar in quantum mechanics, we will 
give the solutions only. There are only bound-state solutions. The energy levels are 

Eq = m n , n r = (33a) 
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£+ ± = ±yfml + 4n r Pfi a , TV = 1, 2, ... , (33b) 

where n r is a radial quantum number. Note that the superscript + for £ indicates the first 
kind of solutions, while the subscript ± indicates the sign of the energy. As before, we have 
negative- as well as positive-energy levels. The corresponding radial wave functions are 



for n r = 0, and 



K r i±\ r ) 



At l±P l e-^L l Z 1/2 (p 2 ) 



.,+ 

n r /± 



for n r 7^ 0, where 



(34) 

(35a) 
(35b) 

(36) 



and L l ^/ 2 (p 2 ), etc., are Laguerre polynomials as employed in Sec. Ill but the argument here 
is p 2 . The normalization constants are determined by Eq. (24) and are given by 



A + 



V2(PPn 



r(/ + 3/2) 



n r 



n r l 



'£: 



r± 



ffln 



C n r ± 



n r 



1,2,.... 



(37a) 



(37b) 



_F(n r + 1 + 3/2) 

It is remarkable that all the above energy levels are independent of the quantum number /, 
and thus all of them are infinitely degenerate. The critical-energy solution (34) is another 
realization of the solution (22). 

Now we consider the second kind of solutions (17). It is easy to show that Eq. (18) gives a 
trivial solution when £ = m n . Thus all nontrivial solutions arise from Eq. (19) which is valid 
for £ ^ m n and in the present case becomes 



d 2 v~ 



2dv~ 
+ -- 



dr 2 r dr 



+ 



£' z -m z n - (21 + 3) /3 fin - p 2 pY 



1(1 + 1) 



v~ = 0. 



(38) 



This is very similar to Eq. (32). The only difference lies in the sign of the third term in the 
square bracket. This difference, however, will render the energy levels quite different from 
those obtained above. As before, we only give the results here. The energy levels are 



£ n± = ±\f< + ( AN + 6 )/^n, N = n r + 1 = 0,1,2,..., 



(39) 



where n r — 0, 1, 2, ... is a radial quantum number and iV is a principal quantum number. The 
superscript — of £ indicates the second kind of solutions, while the subscript ± indicates the 
sign of the energy. The spectrum obtained here has no overlap with that in Eq. (33). The 
corresponding radial wave functions are 



u 



n r Z± ' 



-A- rl± ^^p l+ ^l 2 L^ 2 (p 2 ), 



'N± 



(40a) 
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v~ rl± (r) = A; rl±P l e-^L l + 1/2 ( P 2 ), (40b) 

where p is given by Eq. (36), and n r — 0, 1, 2, . . . is the radial quantum number. The normal- 
ization constants are given by 



A nrl± = (/3/i n )t 



i ,i 



r(n r + I + 3/2) 



n r = 0,l,2,.... (41) 



The energy levels and depend only on the principal quantum number N. Given N, 
I may vary from to N. For a given I, there are 2/ + 2 different solutions. Therefore the 
degeneracy of the level or is 

N 

d N = J2(2l + 2) = (N + l)(N + 2). (42) 

1=0 

In conclusion, in the radially linear electric field, we have two sets of bound-state energy 
levels. The first set is given in Eq. (33), corresponding to the first kind of solutions. The 
second set is given in Eq. (39), corresponding to the second kind of solutions. There is no 
scattering solution here. In contrast, the radially constant electric field studied in Sec. Ill 
admits both scattering and bound-state solutions, though there exists no bound state of the 
second kind. Finally we estimate the "Bohr radius" of the neutron in the present case. This 
is roughly equal to (/3/x n ) _ 2 ; or (3/47rp c /x n ) 2 where p c is the electric charge density producing 
the field. In the MKS system this reads 



3h 



{^4np p c p D J 

where po is the permeability of the vacuum. We take p c = e/a^ where e is the electron charge 
and a is the Bohr radius of the hydrogen atom. For neutrons the above "Bohr radius" has 
the value 4.4 x 10~ 8 m. This is rather small. However, it may be difficult to achieve the above 
electric charge density. 



V. SUMMARY AND DISCUSSIONS 

In the preceding sections we have studied the Dirac-Pauli equation of a neutral fermion 
with nonminimal coupling to a central electric field. By separation of variables in spherical 
coordinates, the stationary Dirac-Pauli equation which involves four partial differential equa- 
tions can be reduced to a system of ODE which involves two coupled first-order ODE for two 
radial wave functions. There are two different kinds of solutions, and thus two independent 
systems of ODE. Bound states of critical energy values can be obtained analytically for a quite 
general class of electric fields, where the degeneracy of the critical energy level turns out to be 
numerably infinite. This reveals the possibility of condensing infinitely many fermions into a 
single energy level. We also discussed a special form of the electric field that supports a finite 
number of critical bound states. Two specific electric fields, one radially constant and the other 
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radially linear, are studied in detail and all the bound-state solutions are obtained in closed 
forms. In the first case bound states exist only for the first kind of solutions, while scattering 
states exist for both kinds. Scattering states are not discussed in detail. In the second case, we 
have two sets of discrete energy levels corresponding to the two kinds of solutions. There is no 
scattering state. It turns out that the energy levels in the first set are all infinitely degenerate. 
In both fields we have negative as well as positive energy levels. Critical energy levels are 
also admitted in both cases, which may be positive or negative depending on the signs of fi n 
and the electric fields. Note that the two critical energy levels are not admitted at the same 
time, however. This spectral asymmetry may likely lead to vacuum polarization after second 
quantization. 

In Sec. II we have shown that the total angular momentum J is a conserved quantity in the 
simultaneous presence of a central magnetic field and a central electric field. But we have not 
discussed the solutions of the Dirac-Pauli equation in this case. In the Dirac representation, 
the stationary Dirac-Pauli equation (6) takes the form 

<r • (p — iji n Ee r )ip = {£ + m n — ji n Bcr • e r )x, (43a) 

<j ■ (p + i[i n Ee r )x = {£ -m n + fi n Bcr ■ e r )(p. (43b) 

These equations are similar to Eq. (8) but more complicated. They are still separable in 
spherical coordinates. We set 

<p(r, 6, 0) = M + (r)/+ (0, 0) + «»/,-(*, 0), (44a) 

X (r, 6, 0) = iv + {r)f hn {0, 0) + *r(r)/£(0, 0)- (44b) 

Substituting these ansatz into Eq. (43) and using the relations (10-12) we obtain the following 
system of ODE for the four radial wave functions. 

-(£ + m n )v + + fi n Bv~, (45a) 

= (S - m n )u + + fJ, a Bu~, (45b) 
= -{£ + m n )v~ + n n Bv + , (45c) 

{£ - m n )u~ + ii n Bu + . (45d) 

If B(r) = 0, one may set u~ = v ~ = which reduces Eq. (45) to Eq. (14) for the first kind 
of solutions, or set u + = v + = which reduces Eq. (45) to Eq. (18) for the second kind of 
solutions. This is what we have done before for a pure electric field. When a magnetic field 
is present at the same time, this is not allowed, however. The essential reason is that K is no 
longer a conserved quantity in this case. All the four ODE in Eq. (45) are coupled to each 
other. It seems difficult to solve them even for a pure magnetic field. We will not go into 
further details of these equations here. 

Let us briefly discuss the nonrelativistic limit of the Dirac-Pauli equation. Consider the 
stationary case with a pure electric field. We can solve Eq. (8a) for x, an d substitute it into 
Eq. (8b) to obtain an equation for ip: 



du + , I , 

— h Unhjv/ v/ = 

dr r 

dv + / + 2 H 
jj, n hv H v 



dr 
du~ 



dr 



+ /j, n hu H u 



dv _ I 

— fX n Ev V 

dr r 
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[a- ■ (p + ifi n E)][cr ■ (p - ^ n E)]y? = (£ 2 - m*)p. 



(46) 



This holds for any value of £ except 8 = —m n , and is valid for noncentral electric field as well. 
To discuss the nonrelativistic limit we consider only positive £ and set 

£ = m n + £'. 

When £' <C m D we get the nonrelativistic limit of Eq. (46): 

[cr • (p + i/i n E)] [cr • (p - i// n E)]</? = 2m n £'<p. (47) 

This has essentially the same form as Eq. (46), and thus the same solutions. However, it 
should be remarked that even when |/x n E | <C m n , Eq. (47) is not valid for those £' comparable 
with m n . For example, in the radially constant field with l/in-^ol <C m n , Eq. (47) is good for 
all bound states, but not for scattering ones with large £, say, £ = 2m n . On the other hand, 
even if |E| is unbounded, Eq. (47) is still valid for small £'. For example, in the radially linear 
field, Eq. (47) may be good for lower levels if \P/J, n \ <C m n . Since Eq. (47) is not simpler, it is 
more convenient to deal with Eq. (46) directly. The nonrelativistic limit with both magnetic 
and electric fields can be similarly discussed, though the situation is more complicated. We 
will not give further details here. 

We have pointed out in Sec. Ill that the radially constant electric field admit scattering 
solutions of both kinds. Though Eqs. (25) and (30) can be solved to give partial wave 
solutions, the scattering problem is difficult to handle in this case since these equations involve 
long-range "Coulomb potentials". An easier situation for the scattering problem may be the 
field E(r) oc r _1 . This will be studied subsequently. 

In this paper we have dealt with (3+l)-dimensional problems. The Dirac-Pauli equation (1) 
has a much simpler form in a (2+l)-dimensional space-time. Indeed, the situation for the AC 
effect is equivalent to a (2+l)-dimensional problem because of the specific field configuration. 
Recently, we have calculated the probability of neutral particle-antiparticle pair creation in 
the vacuum by external electromagnetic fields in 2+1 dimensions, based on the nonminimal 
coupling [14]. Both scattering and bound-state problems in external fields are easier in 2+1 
dimensions. These and other consequences of the nonminimal coupling will also be studied 
subsequently. 
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